Novel type of orbital ordering: complex orbitals in doped Mott insulators 
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An orbital ordering, often observed in Mott insulators with orbital degeneracy, is usually supposed 
to disappear with doping, e.g. in the ferromagnetic metallic phase of manganites. We propose that 
the orbital ordering of a novel type may exist in such situation: there may occur ferro. orbital 
ordering of complex orbitals (linear superposition of basic orbitals d^2_y2 and d,,2 with complex 
coefficients). Despite the perfect orbital ordering, such state still retains cubic symmetry and thus 
would not induce any structural distortion. This novel state can resolve many problems in the 
physics of CMR manganites and can also exist in other doped Mott insulators with Jahn- Teller 
ions. 

PACS numbers: 75.30V, 71.70E 

Mott insulators often display various types of ordering connected with localized electrons. In particular, in the 
presence of an orbital degeneracy, very often an orbital ordering occurs besides the magnetic (spin) one. Classical 
examples are the systems with doubly-degenerate Cg-electrons, e.g. KCuFa or MnFs [1]. At present special attention 
is paid to the manganites Lai_a;A/a;Mn03 (M = Ca, Sr) with the colossal magnetoresistance (CMR). The undoped 
material LaMnOs is a Mott insulator with the orbital ordering. What becomes of this ordering in doped materials 
and which role the orbital degrees of freedom do play in the CMR phase, is a matter of active investigation and of 
hot debate. On the one hand it is usually assumed that the orbital ordering (at least the long-range one) disappears 
in a ferromagnetic metallic (FM) phase (0.16 < x < 0.5) displaying CMR. This is in particular concluded from 
the disappearance of lattice distortion typical for the orbitally ordered states. On the other hand the attempts to 
explain transport properties in this regime using only electron-spin scattering in the framework of the standard 
double-exchange model failed [2] , which forced Millis and coworker to suggest that the orbital degrees of freedom are 
still important in this regime, e.g. producing Jahn- Teller polarons. Nevertheless the orbitals were still assumed to be 
disordered in the FM phase below Tc- 

Until now in all the numerous studies of the cooperative Jahn- Teller effect and orbital ordering only one class of 
solutions was considered: occupied orbitals were always assumed to be certain linear combinations of the basic orbitals 
~ -^(2z^ — — y^) and d^2_y2 ^ '^(^^ ~ v"^) of the form 

\e)^cos^-\z^)+siJ-\x'-y'), (1) 

i.e. linear superposition with real coefficients. Such orbitals give an asymmetric (quadrupolar) distortion of the 
electron density at a given site. Consequently, any ordering involving orbitals of the type (1), be it of "ferro" or 
"antiferro" type, corresponds to a decrease of the symmetry of charge distribution and is accompanied (or driven) by 
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the corresponding lattice distortion. 

There exists however yet another possibihty: there may exist orbital ordering without any lattice distortion. Prom 
the basis set of doubly-degenerate orbitals we can also form a linear superposition of 

the type (1) but with complex coefficients. The simplest example is the normalized functions 



We can form an orbital ordering with the occupied orbitals of such type, e.g. the ferro. orbital ordering (FOO) 
with the orbitals 1+) at each site. It can be easily checked that the electron charge distribution in this state is 
perfectly symmetric, it is the same along x, y and ^;-directions, so that the net symmetry of a crystal would remain 
cubic. Correspondingly, the effective electron hopping matrix elements are isotropic, t^ = t^ = t^ = t* , as well as 
superexchange interaction = jy = = 2t*'^/U. (Here we have in mind the traditional description of the strongly 
interacting electrons by the Hubbard-type model which in case of orbital 

degeneracy is actually a degenerate Hubbard model [1] with the on-site repulsion U and electron hopping between 
different orbitals, to be specified below.) We suggest that just such an ordering is reahzed in the FM phase of CMR 
manganites. If true, this would imply that both the spin and orbital degeneracies are lifted in the ground state in these 
systems which consequently are ferromagnetic in both the spin and orbital variables (but with "strange" orbitals (2)). 

Why have such states not been considered before? Apparently it is connected with two factors. On one hand, in most 
real materials with Jahn- Teller ions studied up to now there always occurred orbital ordering accompanied by lattice 
distortion. On the other hand, in the theoretical description of these phenomena only the operators corresponding 
to an ordering of conventional type appear in the cfFcctivc Hamiltonian (sec below), so that it seemed that other 
possibilities are never realized in practice [3,4]. Most probably, indeed only the "real" solutions (1) can be realized 

in stoichiometric compounds with Jahn- Teller ions. The situation however may be different in doped materials — 
and that is what we propose here. 

It is well known that the electron (or hole) motion in systems with strongly correlated electrons is hindered by an 
antiferromagnetic background and is much facilitated if the background ordering becomes ferromagnetic (see e.g. [5]): 
in nondegenerate Hubbard model it gives rise to Nagaoka's ferromagnetism [6], and in systems with two types of 
electrons — to ferromagnetism due to double exchange. The energy gain in such a ferro. state as compared to the 
antiferro. or paramagnetic (disordered) one is of the order of the electron bandwidth, or of the hopping matrix element 
t per doped charge carrier: AE = —ctx (c is some constant of order 1), whereas the energy loss is of the order of 




(2) 
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J (~ e/U) ■ (1 - x). 

Thus, as argued already in [6], the saturated ferromagnetic state is reahzed for a; > Xc ~ J/t t/U. For smaller x 
the inhomogeneous phase separated state with ferromagnetic droplets in an antiferomagnetic matrix can be realized [7] . 

In exact analogy to this case we should expect that due to the same mechanism a ferro. orbital ordering will be 
established in doped strongly interacting Mott insulators with orbital degeneracy. Which particular orbitals would 
be stabilized at that, may depend on the particular situation; we suggest that it may be the complex orbitals of the 
type (2). 

A convenient way to describe the orbital ordering of double-degenerate Cg-orbitals is to introduce pseudospin 
variables such that e.g. the orbital corresponds to = and the orthogonal orbital (i^2_y2 — to = —5. 
The orbital states (1) considered until now are parametrized by the angle 6 in the (r^, T^)-plane. Cubic symmetry 
is reflected in the ^-plane in ^ symmetry: the state 6 = ^ corresponds to the orbital \x^) ~ (2a;^ — — z'^) and 
= — ^ — to ~ {2y'^ — x^ — z^); these are more or less the orbitals occupied in two sublattices in the undoped 
LaMnOs. When we dope the system with orbital ordering, the motion of the charge carriers would initially (for 
small x) occur on a background of this orbital ordering. Similar to the hole motion on an antiferromagnetic spin 
background, the "antiferro" orbital ordering would hinder the motion of a hole and would reduce its bandwidth. One 
can indeed check that, by making the orbital order ferromagnetic, e.g. occupying the same orbital at each site, for 
instance \z'^) or \z'^ — y^), we would increase the bandwidth and correspondingly decrease kinetic energy of holes; the 
mathematical details are presented below. However such a FOO seems to contradict experimental observations. Thus 
e.g. |^^)-ferromagnetism would lead to a tetragonal distortion of the lattice with c/a > 1 and to a strong anisotropy 
of transport properties; but nothing like that is observed experimentally. Probably that is why this possibility, which 
by analogy with the spin case seems quite natural, is never considered for CMR phase of manganites, and it is usually 
assumed that the orbital ordering 

is simply lost in the FM phase leading e.g. to an orbital liquid [8] However if at such a FOO the complex orbitals (2) 
would be occupied — there would be no contradiction with the structural data, and still the motion of doped charge 
carriers would be unhindered. This is the main idea of the present paper. In the pseudospin language introduced 
above the state (2) is an eigenstate of the operator which also exists in the algebra of r = | operators. 

For strongly interacting electrons with one electron per site one can describe the effective spin and orbital interaction 
by the Hamiltonian which schematically has the form [1] 
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W = ^ Js{SiSj) + JrinTj) + JsT{SiSj){TiTj) , (3) 

where the first and the third terms are due to superexchange, Jexch ^ Js ~ Jst f^/U, and in the second term also 
the Jahn- Teller interaction of degenerate orbitals with the lattice contribute, Jr = Jexch + Jjt The interaction (3) 
is in general anisotropic with respect to r-operators and usually it docs not contain terms with t"^ . Consequently 
r^-states do not appear in the mean-field approximation which is nearly always used to treat insulators with orbital 
degeneracy. But it does not mean that such states cannot appear in certain situations, e.g. for doped systems. 

Using the specific form of the Cg-orbitals and of the corresponding hopping integrals between different orbitals 
in different directions, one can easily calculate the one-electron band structure for different types of orbital order- 
ing. Denoting the hopping between ^;^-orbitals for a pair along ^;-direction as t, we have: t^2 ^2 = t, t^Jz^^ = i/4, 
= 3i/4, = ±V3t/4, i^,^^, = ejy2 = t/2, etc., see e.g. [1] 

Generally speaking, in contrast to spin case (see [5] for discussion and references), an antiferro. orbital ordering, 
e.g. that in LaMnOs, does not completely suppress the hole motion even if one ignores quantum effects: the hole 
can always hop to a neighbouring site without destroying orbital order along its trajectory. This is connected with 
the fact that the nondiagonal hopping in orbital channel is in general allowed, and pseudospin projection is not 
conserved during hopping: the nondiagonal hopping integrals tz2 ;^2_y2 are in general nonzero. Thus the hole can 
here move without leaving a "trace" of wrong spins [9] and, consequently, there will be no "confinement" . However 
the bandwidth of this coherent motion without disturbing the background orbital ordering is reduced as compared to 
the bandwidth for FOO. Thus, using the values of hopping integrals t^^-'J given above, one obtains for the undoped 
LaMnOa (alternation of x"^ and orbitals) the spectrum 

e{k) = —2 [if(cos kx + cos ky) + \t cos kz\ (4) 

so that the minimum energy of the hole (the bottom of the band) will be Emm = £(fc=0) = — 2.5t. If however we 
would make a FOO, e.g. occupying at each site z^-orbital, the spectrum will be 

£^^-/erTO(^) ^ _2t ^1 (cos kx + COS ky) + COS k^] (5) 

2 J- 

and s^jnln" — Thus we indeed see that the energy gain in the FOO state is of order ^ t. The same minimum 
energy is reached also for [x^ — y^)-ferro. ordering and for FOO with any orbital of the type (1). From (1) one can 
easily obtain that 
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i|,e = (^l*'l^)=*cos2|; t^// = |(cos|±x/3sin|)2 , (6) 

And one obtains that the bottom of the spectrum with the hopping integrals given by (6) does not depend on 9 and 
coincides with the value —3t. 

But exactly the same Smin is also reached for the r^-ferromagnctism, where the state (2) is occupied at each site. 
Using the values of t's presented above, one obtains that try,TV = \t in all three directions, and 

^T^-ferrof^j^^ = -t(cOS + COS ky + COS k^) (7) 

SO that the spectrum is indeed isotropic, and the minimum energy is also equal to —M. 

Note that, similarly to the nondegenerate Hubbard model, we can use the simple dispersion relations (4), (5) only 
for the electron motion which does not destroy the background ordering. The energy gain in the FOO state is again, 
similarly to the Nagaoka case [6], ~ tx, and the energy loss is ~ J (in orbital sector ~ + Jst)- Here again we obtain 
the totally ferromagnetic state both in spins and in orbitals for x> Xc'^ J/t, and one can have phase separation for 

As we saw above, for small x different FOO states are equivalent as to the gain in kinetic energy. On the other hand 
we can see that r^-ferro. state is favourable as to the loss of the effective exchange energy (3). Indeed, the effective 
Hamiltonian (3) contains the orbital operators in combinations rf r? , rf rj and T^Tj with positive (antiferromagnetic) 
coeficients J^, Jst [!]• In the ground state of the undoped system a certain spin and orbital order is realised which 
minimizes the total energy — and it will be an antiferro. ordering of real orbitals (1). If however we force our system 
to be ferromagnetic both in spins and in orbitals, e.g., as argued above, by doping, we strongly increase this exchange 
energy — and we can minimize this energy loss by chosing FOO of r^-type. 

One can check directly that the superexchange part of the energy is indeed smaller in the r^-FOO state as compared 
e.g. with the FOO state of real orbitals. Thus, taking the standard expression for the exchange integrals J = 2t^/{7 
and calculating exchange constants with the proper values of t, we obtain e.g. that for the z^-FOO the energy of the 
ferromagnetic state (per site) in mean-field approximation can be written in the form: 

p.2^-FOO _ ^ I r" (a.\ 

^ferro ~ 4 [/ ' ^ ' 

where C = —3^. The same value of Efg^ro we would obtain for (a;^ — y^)-FOO and for all other FOO states with 
arbitrary state |^) of the type (1): again one can easily show that with the values of hopping integrals tgg (6) the 
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magnetic energy of both ferro- and antiferromagnetic states does not depend on 9 and is given by (8). Similar 
treatment of the r^-FOO state shows however that the net exchange energy in this state is lower: 

^ferro ~ 2U 

where C is the same as in Eq.(8). Similarly, the intersite Jahn- Teller energy in FOO state with real orbitals is positive, 
^ JjT, whereas it is zero for the r^-FOO state; again from this point of view the r^-FOO state wins against 
all other possible FOO states with real orbitals. 

As the electronic energy of doped charge carriers coincides for the states (1) and (2), and the magnetic energy 
of the latter is lower, one can conclude that the proposed ferromagnetic state with r^-FOO (occupation of complex 
orbitals (2)) will be the best among all possible ferro. ordered states — and these, we argue, may be stabilized 
by the usual Nagaoka mechanism. This state can possibly be realized in the ferromagnetic metallic phase of OMR 
manganites. The finite band filling would increase electron kinetic energy somewhat faster for r^-ordering as compared 
e.g. to z^-FOO or (x^ — i/^)-FOO (the energy spectrum in the latter cases is anisotropic and the density of states at 
the band edge is higher, cf. [10]). Consequently the r^-FOO may be destabilized at still higher doping levels. 

The state with r^-FOO need not be necessarily homogeneous: as in the pure spin case, for small doping there may 
exist a tendency towards phase separation due to an instability of the homogeneous canted spin state [7,11]. But, as 
follows from the arguments given above, the structure of the hole-rich regions can well be of r^-FOO type. 

^From the point of view of symmetry the proposed r^-FOO corresponds to the A2g (or r2)-representation of the 
cubic group [12]. The corresponding order parameter is 

Txyz = {Slxlylz) ) (10) 

where ( ) is the ground state average, <S denotes 

symmetrization and Ix, ly, Iz are the corresponding components of the momentum operator / — 2. This follows 
from the group symmetry, and we checked by direct calculation that it is the lowest nonzero average in our state. 
One can show that in the basis \z^), \x'^ — xp') T^yz is indeed proportional to the r^-Pauli matrix. From eq. (10) we 
see that the order parameter in our state is a magnetic octopole, cf. [3,4]. It is in principle possible, although not 
easy, to observe such ordering experimentally. Probably the most promising would be the resonance experiments like 
NMR or circular X-Ray dichroism. One can also show that this state should have a piezomagnetism. 

It is not apriori clear whether this r^-FOO and the ferromagnetic spin ordering would occur at the same temperature. 
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The direct product of corresponding irreducible representations A2 and Ti does not contain unit representation Ai 
[12], so that there will be in general no linear coupling between these order parameters. However if the rhombohedral 
distortion T2g, present in manganites at high temperatures (S. W. Cheong, private communication), would still exist 
at low temperatures, then there will exist linear coupling of these three types of ordering {A2 x Ti x T2 e j4i) and the 
r^-FOO and the ferromagnetic ordering will occur at the same Tg. This could help to resolve the problem pointed 
out in Ref. [2], that the change of resistivity through Tg in purely spin case is too small: one would get an extra 
scattering in a "para" phase — an orbital-disorder scattering — in addition to the spin-disorder one. (Alternatively 
one can speak of the increase of resistivity above not due to orbital fluctuations, but, again in analogy with the 
spin case [9,13], due to the effective narrowing of the band in the orbitally disordered phase.) 

Summarizing, we suggest in this paper that, similarly to magnetic ordering, doping of orbitally degenerate Mott 
insulators containing Jahn- Teller ions could stabilize the ferro. orbital ordering of special type, the occupied orbitals 
being "complex" — i.e. the linear superpositions of basic orbitals dx2_y2 and (1^2 with complex coefficients. Such 
orbitals, despite perfect ordering, do not induce any structural distortion. At the same time the motion of charge 
carriers on such ordered background is completely free, and there is no extra band narrowing; it is just this factor 
that favours such ferro. orbital ordering. Thus the origin of the ferro. orbital ordering in doped systems is the same 
as that of the Nagaoka ferromagnetism in a partially-filled Hubbard model. The ferro. ordering of complex orbitals 
gives the same band energy as the ordering of conventional real ones but has lower exchange and Jahn- Teller energy; 
therefore the ferro. ordering of complex orbitals may be preferable. This idea may explain the main properties of 
the colossal magnetoresistance manganites in the most interesting ferromagnetic metallic concentration range: in this 
picture they are perfectly ordered with respect to both spin and orbitals, but with the complex orbitals occupied at 
each site. In particular, this could help to explain the sharp drop of resistivity below T,,, with the crystal structure 
remaining undistorted. Similar states may in principle exist also in other doped Mott insulators with Jahn- Teller ions 
which constitute a large interesting class of magnetic materials with very specific properties. 
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